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We investigate the ensemble inequivalence of the spherical spin glass model with nonlinear 
interactions of polynomial order p. This model is solved exactly for arbitrary p and is shown 
to have first-order phase transitions between the paramagnetic and spin glass or ferromagnetic 
phases for p > 5. In the parameter region around the first-order transitions, the solutions give 
different results depending on the ensemble used for the analysis. In particular, we observe 
that the microcanonical specific heat can be negative and the phase may not be uniquely 
determined by the temperature. 
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1. Introduction 

Statistical ensembles play key roles in the investigation of systems with a large number 
of particles. The most fundamental one is the microcanonical ensemble, and it is believed 
that we can derive other ensembles from the microcanonical ensemble. !) Alternatively, we 
can relate the microcanonical ensemble to the canonical ensemble by the Legendre transfor- 
mation. !) It implies that the physical properties of a system do not depend on the statistical 
ensemble in the thermodynamic limit. In other words, these ensembles are generally consid- 
ered to be equivalent. This fact enables us to choose a statistical ensemble which is useful for 
practical calculations, and the canonical ensemble is usually the first choice. However, this 
property is known to hold only for short-range interacting systems. For long-range interacting 
systems, we have to reconsider whether or not ensembles are equivalent. 

In fact, the phenomena that the physical properties depend on statistical ensemble emerge 
in systems with long-range interactions. In particular, for the systems with first-order tran- 
sitions, ensemble in-equivalence has been reported in previous studies. 2 ~ 16) An interaction 
decaying in a power low r a (a < d), where r is the distance between particles and d is the 
spatial dimension, is said to be long ranged. 5 ~ 8) Long-range interacting systems with first- 
order phase transitions have a striking feature that the specific heat can be negative under 



the micro-canonical ensemble. Such ensemble inequivalence has traditionally been discussed 
mainly in astrophysics. 9 ~ 12) 

Recently, these problems outside of astrophysics have been analyzed from the view 
point of ensemble inequivalence, results of which have further promoted the understand- 
ing of ensemble inequivalence. Many non-trivial results related to ensemble inequivalence 
have been reported in spin systems without disorder 13_17) and driven systems with local 
dynamics. 18, 19) Analyses of long-range interacting spin systems with disorder have only re- 
cently been initiated, 20 ~ 23) but our knowledge on ensemble inequivalence in systems with spin 
glass transitions is still limited. 20 ' 22 23) Therefore, it is an interesting problem whether or not 
ensemble inequivalence are observed in the other spin systems with disorder. 

In this paper, we analyze the spherical spin glass model with non-linear interactions both 
by the canonical and microcanonical ensembles, to investigate whether or not ensemble in- 
equivalence emerge in this spin glass system. This model is a generalization of the linear 
model introduced by Kosterlitz et al 24 ' 25) in conjunction with the nonlinear interaction for the 
non-disordered cases 17 ' 26 ~ 29 ) and has the advantage that all calculations can be done exactly. 
We can thus avoid the replica method and do not have to worry about the limit of applicability 
of the replica method. 

The organization of this paper is as follows. The model is defined in § 2. In § 3, we 
calculate the partition function in the canonical ensemble (the number of states in the micro- 
canonical ensemble) and obtain the free energy density (the entropy density). We analyze the 
behavior of thermodynamic functions and phase diagrams in § 4. The last section is devoted 
to summary and conclusion. 

2. The Model 

Let us introduce the spherical model of spin glasses with non-linear interactions. It is 
defined by the Hamiltonian 




(1) 



Here, V (x) is given by 



V(x) = 



1 



[a + *y-i], 



(2) 



2P~ l p 



and gij is a dimensionless random parameter, whose probability distribution P (^J is Gaus- 
sian with average g /N and variance 1 /N, 




The spin variables S , obey the spherical constraint 

N 

S 2 = J]s 2 =N. (4) 

1=1 

In the case of p = 1, V(x) is equal to x and the model reduces to the spherical model 
introduced by Kosterlitz et al. 24) 

In the absence of disorder, 26-250 the system shows first-order transitions and ensemble 
equivalence for p > 5 as long as the interactions are short-ranged. 17) 

3. Free energy and entropy 

In this section, we calculate the partition function and the number of states, from which 
we derive the canonical free energy and the microcanonical entropy. 



3.1 Free energy in the canonical ensemble 

First we solve the model in the canonical ensemble. The partition function of the system 
can be written as 



Z = Tr exp 



NfiV 
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s(s 2 -n), 
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where {3 is the inverse temperature and the trace denotes integrations over the spin variables. 
We introduce the auxiliary variables z, p and A to facilitate the integration, 



dzdpdA Tr exp 



N 



N 



. (=i 



-A 



dzdpdA exp 



N/3V (p) + Nz- NAp + log Tr exp 



Np-^ijSiSj 

N 

- Z J]S 2 + A^ U S { S j 



i=i 



KJ 



(6) 



By diagonalizing the random matrix we calculate the trace part as a Gaussian integral. If 
we denote the eigenvalues of the random matrix |^ (; } as {^ a }, we can write 



dzdpdA exp < iV 



fiV(p) + z-Ap-^J]\og(z-^^ 



(7) 



The saddle point equations for the evaluation of the integral for large N are 

dV 

* = fi—, (9) 
dp 

In the thermodynamic limit, the summation over a is replaced by an integral, 

i r°° 

kJ]^ dfa/Jiifa). (11) 

The density of eigenvalues //(£ a ) can be expressed as 

Ho(£ a ) (ft < 1) 

= i . (12) 

Mta) + J (ft > !)• 

Here /i (ft) obeys the semicircular law a/4 - £l/2n. 30) If ft > 1, the eigenvalue spectrum is 
modified by an isolated eigenvalue £ m = £ + 1 /ft- 31) The thermodynamic function for f > 1 
describes the ferromagnetic (FM) phase. 

Using eqs. (11) and (12), we evaluate the saddle point equations, eqs. (8) and (10), to 
determine the state of the system. We eliminate the auxiliary variables A and p from the saddle 
point equations. Then the free energy density f = -T log Z/N and the inverse temperature f3 
of each phase are given as functions of z, from which / is determined as a function of /?. Note 
that we should carefully treat the estimation of saddle point equations when we replace the 
summation of eigenvalues by an integral. Details are described in Appendix. 

The paramagnetic (PM) and spin glass (SG) phases exist when the ferromagnetic bias is 
small, £ < 1 . The free energy density and the inverse temperature of the PM phase are 



-/Wpm = fcrnV^—^] - \z + l - log2 + 1, (13) 



(14) 



(l +y/2z~=T/2f 

which are valid for 1/2 < z < 1. It is seen in eq. (14) that z = 1/2 corresponds to {3 = 
(T —* oo), which justifies us to identify this range of z with the PM phase. The SG phase 
exists when z > 1 . The free energy and the inverse temperature are given by 

-£sg/sg = ^sGVll-^)-ilogz+ilog2 + i (15) 



(16) 



(2-l/2zy- r 

We can confirm that the free energy densities of the PM and SG phases take the same value 
at z = 1, which agrees with the previous result for p = I. 24 ' 25) 

The free energy density and the inverse temperature of the FM phase for £ > 1 are given 

by 
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Notice that the free energy density and the inverse temperature of the FM phase are equal to 

those of the SG phase when £ = 1 corresponding to £ m = 2. From this fact, we confirm that 

/ (T) is a smooth, continuous function at transition point between the FM and SG phases for 

arbitrary p. 



3.2 Entropy in the microcanonical ensemble 

We next calculate the number of states for the microcanonical ensemble as described in 
Appendix, 

Q = Tr6(E-H) 



= dtdzdpdA exp In it(e + V (p)) + z - Ap - log {z - ^A£ a J 



,(20) 



where e is the energy density E/N. The saddle point equations of the microcanonical ensem- 
ble have very similar expressions to the canonical case in eqs. (8)-(10), 

dV 

A = it—, (22) 
dp 

P = 2^2^' (23) 
with an additional equation derived by the derivative of the exponent of eq. (20) with respect 
to t, 

= e + V(p). (24) 



We solve these equations and find the microcanonical entropy density s = log Q./N as de- 
scribed in Appendix. The PM and SG phases are defined for £ < 1- The entropy and energy 
densities of the PM phase for 1 /2 < z < 1 are given by 

1 1 3 
spm = --z+-log2+-, (25) 

*m = -vi^El). (26) 
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The entropy and energy densities for z > 1 describe the SG phase, 

1 1 1 

*sg = --logz+-log2+-, (27) 

esc = - y ( ! -^)- ( 2g ) 

Similarly, the entropy and energy for the FM phase for £ > 1 are given by 

1 1 1 

Sfm = ~^ogz + -log 2 + - +h(t m ), (29) 



2 1 2z 



(30) 



Here h (g m ) is the same function as in eq. (19). The FM solutions are equal to the SG ones at 
^o = l- Therefore, s (e) continuously changes between the FM and SG phases. 

4. Comparison of the thermodynamic functions and the phase diagrams 

In this section, we compare the canonical and microcanonical solutions for p = 1 to 
p = 5. We find that the phase transition between the PM and SG phases is of first order in the 
canonical solution for p > 5, which results in ensemble inequivalence. 

4. 1 Thermodynamic functions 

First, we focus our attention to the PM and SG phases (£ < 1)- In Fig. 1, we plot the 
inverse temperature /? PM and /? S g of eqs. (14) and (16) for the canonical ensemble in Fig. 1(a) 
and the energy density e PM and 6 S g of eqs. (26) and (28) for the microcanonical ensemble 
in Fig. 1(b) for p = 1, 2, 3, 4 and 5. Notice that /3 PM is for z < 1 and /3 S g is for z > 1 in 
Fig. 1(a) and similarly for Fig. 1(b). In the canonical case,/?(z) is a monotonically increasing 
function of z for p < 4. When p = 5, in contrast, /? (z) is not monotonic and z is not uniquely 
determined for a given f5 in a certain range. We can understand this fact as an indication 
that the transition between the PM and SG phases is of first order for p = 5. We also find 
that the inverse temperature of p > 6 behaves qualitatively in the same way as the one for 
p = 5, which we do not show in the figure for simplicity. By contrast, e (z) is a monotonically 
decreasing function of z for any p as seen in Fig. 1(b). We also observe that first-order phase 
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(a) The canonical inverse temperature (b) The microcanonical energy density 

Fig. 1. The inverse temperature in the canonical ensemble in Fig. 1(a) and the energy density in the micro- 
canonical ensemble in Fig. 1(b) for p = 1, 2, 3, 4 and 5 as functions of z. When p < 4, ft (z) is a monotonically 
increasing function. e(z) is a decreasing function for any p. 



transitions appear between the PM and FM phases for p = 5 and £ > 1, because ft (z) is not 
a monotonically increasing function of z (Fig. 2). The same is true for p > 5 but not p < 4. 
As can be understood from these results, in order to discuss the ensemble inequivalence, it is 
sufficient to pick up the typical cases of p = 1 and p = 5. 
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Fig. 2. The inverse temperature for p = 5 and £o = 11, 1.3, 1 .5, 1.7 and 1 .9 as a function of z- 




Next we consider the thermodynamic functions in the PM and SG phases (i.e. £b < 1)- 
We plot the free energy density / (T) of eqs. (13) and (15) as a function of T in Figs. 3(a) and 
3(b) for p = 1 and p = 5. For p = 1, the critical temperature is T c = 1 as seen in eqs. (14) 
and (16) with z = 1. We confirm that the phase transition is of second order because d 2 f/dT 2 
is discontinuous at z = 1. For p = 5, as we can see in Fig. 3(b), the free energy density f(T) 
shows a first-order phase transition at T c = 0.31, which has already been expected in the 
previous discussion on /3 (z). 
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(a) The microcanonical inverse temperature for p = 1 (b) The microcanonical entropy density for p = 1 



Fig. 4. The inverse temperature /3 (e) (Fig. 4(a)) and the entropy density s (e) (Fig. 4(b)) in the microcanonical 
ensemble for p = 1. A second-order phase transition exists between the PM and SG phases at e = -0.5. 



We also study the entropy s (e) and the inverse temperature /3 (e) for p = 1 and p = 5 under 
the microcanonical ensemble. As the inverse temperature f3 (e) for p = 1 is a monotonically 
decreasing function in Fig. 4(a), the shape of the entropy density s(e) is concave (Fig. 4(b)). 
The phase transition between the PM and SG phases is of second order because ds/dT shows 
a singularity at T = 1. In this case, ensemble inequivalence does not appear and we are 
able to obtain the free energy, which is the same one as the canonical case, by a Legendre 
transformation of the microcanonical entropy. On the other hand, when p = 5, the energy 
density e has multiple solutions for a fixed value of jS in a certain range as seen in Fig. 5(a). 
In this region, the entropy becomes a nonconcaved function (Fig. 5(b)) and the specific heat 
takes negative values, which is a clear indication of ensemble inequivalence. In contrast to 
previous studies 21 ~ 23) , the microcanonical entropy for p = 5 smoothly changes from the SG 
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(a) The microcanonical inverse temperature for p = 5 



(b) The microcanonical entropy density for p = 5 



Fig. 5. The inverse temperature /3 (e) (Fig. 5(a)) and the entropy density s (e) (Fig. 5(b)) in the microcanonical 
ensemble for p — 5. j3 (e) is not a monotonically decreasing function and s (e) is not concave. 



solution to the PM solution at e = -0.82. 
4.2 Phase diagram 

Let us next analyze the phase diagram on the (T, £ ) plane for p = 1 and p = 5. For p = 1, 
it is well known that only second-order phase transitions exist. 24 ' 25) In this case, we do not 
observe ensemble inequivalence. 




Fig. 6. The canonical (CE) and microcanonical (MCE) phase diagrams are denoted in the black solid line and 
the red dashed line, respectively, for p — 5. The PM phase exists in the upper part. The SG and FM phases exist 
in the lower part for < 1 an d > 1> respectively. In the microcanonical phase diagram, the phase is not 
uniquely specified for a given value of T around the boundary between the PM and SG/FM phases. 



By contrast, differences emerge for p = 5 as shown in Fig. 6. In the canonical ensemble 
(the black solid lines), the phase transitions between the PM and SG or FM phases is of 
first order as explained § 4.1 and the phase transition between the SG and FM phases is 
of second order. In the microcanonical ensemble (the red dashed lines), there are regions 
where the phase is not determined uniquely by T which is not the control parameter of the 
microcanonical ensemble. For £ < 1, as the energy density is decreased, the temperature 
decreases monotonically down to T = 0.29. When the temperature reaches the value on the 
lower red dashed line (the line at T = 0.29), the temperature turns to increase to T = 0.32 
(the upper dashed line) and the state changes from the PM phase to the SG phase at T = 
0.3 1 . This non-monotonic behavior of temperature means that the specific heat takes negative 
values. After the temperature reaches 0.32, the temperature turns to decrease again down to 
T = 0. Such a result is observed also between the PM and FM phases for £ > 1 • These 
observations mean that the state of the phase is not determined by the temperature in the 
range 0.29 < T < 0.32 on the microcanonical phase diagram for £ < 1, and similarly for 
£ > 1. This behavior has been pointed out in other cases as reported in refs. 21, 22 and 23. 

5. Conclusion 

We have derived the canonical free energy and the microcanonical entropy of the spherical 
spin glass model with nonlinear interactions and have shown that ensembles are not always 
equivalent. In particular, we have studied the behavior of the canonical free energy and the 
microcanonical entropy for p = 1 and p = 5. For p = 5, the microcanonical entropy is not 
concave in a certain parameter range. As a result, the specific heat can take negative values. In 
addition, we analyzed the phase diagrams for p = 1 and p = 5 and have shown differences in 
the phase diagrams on the (T, £ ) plane for p = 5 between the canonical and microcanonical 
ensembles. 

An important advantage of our study is that we have exactly solved the highly non-trivial 
model of spin glasses both by the canonical and the microcanonical ensembles without re- 
course to the replica method. We have thus unambiguously established the presence of en- 
semble inequivalence when quenched disorder exists in long-range exchange interactions. 
We expect that our present work paves the way toward further understanding of ensemble 
inequivalence in systems with disorder. 



Appendix: Derivation of the Free energy and Entropy 

In this appendix, we derive the free energy densities and the inverse temperatures of the 
canonical ensemble, eqs. (13) - (19), from the saddle point equations (8) - (10) and also obtain 
the entropy densities and the energy densities, eqs. (25) - (30). 

First, we derive eqs. (13) and (14). We evaluate the integral of the saddle point equations 
(8) and (10) by using the density of eigenvalues [i (% a ) for £ < 1 and obtain 

l = ^[z-y/z 2 -A 2 ], (A.l) 

9 = 2^ [ 2Z ' ~ 2Z ~ ^ • (A - 2) 
When A < 1, we can solve eq. (A.l) for A and express p as a function of z, 

A = ^J2^T, p = ^2^1, (A3) 

and obtain the inverse temperature J3(z) of eq. (14) by inserting eq. (A. 3) into eq. (9). We 
determine the range of z by the restriction A < 1 as 1/2 < z < 1. In this range, the summation 
on the right hand side of eq. (7) is estimated as follows, 

= z-log2-i (A.4) 

Thus, the free energy density and the inverse temperature of the PM phase are given by eqs. 
(13) and (14). 

Next, we consider the saddle point equations for A > 1. To derive the free energy density 
and the inverse temperature in the SG phase, we should examine the appropriateness of inte- 
gral approximation of eq. (11). As in the case of Bose-Einstein condensation, we separate the 
contributions of the largest eigenvalue i; a = 2 and the rest in eq. (8) as 



1 = 1—!— 



N2z-AZ ao N&Tz-Xa 



_________ r^j^t. (a.5) 

2Nz-A J_ 2 qa 2z-A% a 

Let us suppose that the denominator of the first term in eq. (A.5) is of order 1 /N, 

z-A,±- r -A. (A.6) 

Then the saddle point equation, eq. (8), for A > 1 (i. e. z > 1) is written as 

l=c+ df a — — — =c + - -c + -. (A.7) 
J_ 2 2A-A$ a A z 



In a similar way, we evaluate eq. (10), 



1 £ao 1 V 



2^ - A£ a 

= c+ 2l^-T z (A - 8) 

Here we eliminate c using eq. (A. 7) to express p as a function of z. We obtain the inverse 
temperature in the SG phase of eq. (16) by substituting eqs. (A. 6) and (A. 8) into eq. (9). To 
derive the free energy density in the SG phase, eq. (15), we estimate the summation term in 
the partition function using z = A, 

Jj Yj log [ A ~ \^°) ~* J d ^ a M) ^ log [ A ~ 

^logz-log2 + i, (A.9) 

where we notice that the contribution from the term of £ ao = 2 vanishes in the thermodynamic 
limit. As a result, the free energy density and the inverse temperature in the SG phase are 
given by eqs. (15) and (16). 

Finally, we derive eqs. (17) - (19) by taking account of the ferromagnetic bias. When 
£ o > 1, the isolated eigenvalue % m = £ + 1 /£o appears. The saddle point equation of eq. (8) in 
the FM phase is written as 

1 = <% a - — 

J_oo lZ - A^a 
1 1 



f d% — — 

J 2 a 2z - A% a 



+ (A. 10) 



N2z-A£ n 

The ferromagnetic ordering appears as a result of the first term in the final expression of 
(A. 10). Assuming that 

2z = ^ m + ^^A^ m , (AM) 

we can evaluate eq. (A. 10) as 

1*0 (&) 



f-2 ^"^ m -A^ a 

= d+ f z [^ m ~^~ 4 \- (AA2) 



The saddle point equation of p, eq. (10), for the FM phase is 



i r°° 

fm , 1 f* ii. \ fa 

- J df a yU() (£/ 



1 



^2 



2/V 2z - Ag„ 



fn 



1 



(A.13) 



where we have eliminated J using eq. (A. 12) to obtain the final expression. The inverse tem- 
perature of eq. (18) is given by plugging eqs. (A. 11) and (A.13) in the saddle point equation 
(9). We can obtain the free energy density eq. (17) using 



^logz-log2 + --2fc(£n). 



(A.14) 



Next, to show eq. (20), we write 

Q= J lt. tre i<(E-») d ( N -S 2 ) 

= hf dtdzdp tr exp | zY ( £ + nv &v + Nz - z J] s ^ 6 Np ~^f ^ ijS iS j 



2nJ 



dtdzdpdA exp</V 



it(e+ V(p)) +z-pA - — logtrexp 



-z^ + A^jSiSj 



V i 



(A.15) 

and which is eq. (20) when we diagonalize random matrix and carry out the spin trace 
as a Gaussian integration. 

The entropy densities and the energy densities in the microcanonical ensemble, eqs. (25) 
- (30), are derived similarly as before by the observation that it plays a similar role to /3 in 
the microcanonical case. Notice that we do not have to express it as a function of z in the 
microcanonical case, because the control variable of the system is not it but e. Therefore, we 
only have to obtain e as a function of z using eq. (24) to derive the entropy densities and the 
energy densities of eqs. (25) - (30). 
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